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<^ ■ Abstract 



Generalized symmetry integrability test for discrete equations on the square lattice is 
studied. Integrability conditions are discussed. A method for searching higher symme- 
tries (including non-autonomous ones) for quad graph equations is suggested based on 
characteristic vector fields. 
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^ 1 Introduction 

q 

■^J- \ Consider nonlinear equations on the quad graph (or double discrete chains) of the form 
O 

= f(u,Ui,Ui). (1) 

Here the unknown u = u(m,n) is a function of two discrete variables m,n. For the sake 
of convenience we use the following notations: Uk = u{m + k, n), uu = u(m,n + k), U\ t \ = 
u(m + l,n + 1). Function / is supposed to be locally analytic, it depends essentially on all 
three arguments. In other words equation (TjQ) can be rewritten in any of the following forms 

= f hj (u,Ui,Uj), with i = ±1, j = ±1. (2) 

Equations of the form ([1]) have important applications and the problem of the complete 
description of their integrable cases is challenging. At present there are various approaches for 
studying integrability of discrete models. In [TJ, [2], [3], 0], the property of consistency 
around a cube is considered as a criterion of integrability. Another classification method called 
vanishing of the algebraic entropy is developed in [6], [7], [8], [9]. The symmetry approach is 
approved to provide a powerful classification tool for integrable discrete and continuous models 

no], mi, n2i, ca, nm, m 
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and 



It was observed earlier that characteristic vector fields provide an effective classification 
tool for integrable models either in the case of Darboux integrability [18], [23], [21] or S- 
integrability [19], [20], [21]. Recall that in the discrete case these vector fields are defined as 
follows Yi = D~ x -^r-D and Y_\ = D^—^D~ l and have the following coordinate representation 

d d 1 d did 
Yi = — + x- 1 Vxx x - 1 \- (3) 

OU OUi X-x OU-i OU 2 X„iX_ 2 OU-2 

d d 1 d did 

Y-i = tt + Vj: h +yyi^~ + j; + 4 

du ou\ y-i du-i du 2 y-iy-2 ou-2 

where x = ^-i( gZfeHiM ) = - ffi'^"'" 1 '^^ and y = D( dfl '^ u ' u ^-^ ) = - ff^-.y . 
Here the shift operators D and 5 act due to the rule Dh(m,n) = him, + l,n), Dh(m,n) = 
h(m,n + 1). Vector fields Y±, Y_i, X 1 = -^—^ and X_i = g|- considered as operators applied 
to the variables Ui, u_i, u, u±i, u± 2 , ... satisfy the following conjugation relations 

DX 1 D~ 1 = pXu DY X D- 1 = -Yi, DX^D' 1 = qX_ u DY^D' 1 = -Y_u (5) 

x y 

where v = _p( , d/- 1 -- 1 K»-i^-ih = 1 anc [ Q = jj( Mll^^lMl ) = 1 

" V du—i I df 1 ' 1 {u,u\,u— i)/du— i " ^ dui ' df(u,ui,u—i)/dui' 

To stress the close connection of these operators with the symmetry approach we note that 
restrictions of the characteristic vector fields Y\ and Y-\ on the set of functions depending 
only on coincide with the differential operators, introduced by P.E.Hydon [TU] and 

successfully applied in [14] , |llj . [T2] to look for conservation laws and symmetries 

d f 1 ' 1 d f' 1 ' 1 d 
a _ £j£ Ju (C\\ 

du f^- 1 d Ul fc}r l 9u-i [ ) 



and 



B = - - ^— - l^J- (7) 

du f Ul dui full 1 du -i 



To compare operators (J4]), (jSJ), ([7]) use the following formulas 

l/z-i = -fu 1 '- 1 /^ 1 and l/y- 1 = -f- 1 < 1 /f-_\ 1 . (8) 

Our work was stimulated by the results obtained by Levi and Yamilov. In the over mentioned 
articles [TT], [T2] they made a good progress in a formalization of the theory of integrable discrete 
models and gave an effective tool to look for autonomous symmetries of ([I]) of the form 

u t = g(u, Mi, u-i, Hi, u-x). (9) 

The goal of the present article is to suggest a method suitable for searching higher sym- 
metries of equation ([I]) of general form f[T3"j) depending also on discrete arguments m, n. It is 
important since there are integrable equations ([T]) having no any symmetry of the form ([9]) but 
having symmetries ( !T3|) . We know at least two examples of such a kind, one of them is equation 
( ITTj) considered below and the other is Adler's discretization of the Tzitzeica equation (AdT) 
recently suggested in [25] . 

In [TB] an alternative approach to discrete equations is studied. Necessary conditions of 
integrability are found under assumption that equation ([T]) admits a recursion operator given 
as a formal series of a special form 

R = r x D + r + r^iD' 1 + r_ 2 D~ 2 + ... , 
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with coefficients Tj. 

The article is organized as follows. In section 2 we discuss on a sequence of integrability con- 
ditions implied by the requirement of existence of symmetries of sufficiently high order. These 
integrability conditions are derived by consecutive differentiation of the linearized equation with 
respect to the dynamical variables. They are in fact functional equations. We give the first four 
of the conditions in an explicit form (see Proposition 1 below) . Among them two were deduced 
earlier in [TT]. It makes a sense to compare these integrability conditions with those derived in 
[16J. It is an intriguing question, whether AdT equation satisfies the integrability conditions 
from [16] . In the next section 3 by using characteristic vector fields we deduce differential con- 
sequences of the integrability conditions from section 2. The consequences are given as systems 
of linear first order PDE. Compatibility conditions of these systems are effective necessary in- 
tegrability conditions since they are easily checked for any given quad graph equation of the 
form ([1]). At the same time the systems provide an algorithm of constructing symmetries. In 
Lemma2 it is proved that any symmetry (jTBl below splits down into a sum of two functions like 
g(m,n,Uj,Uj-i, —,Uj>,u k ,u k -i, —u k >) = F(m,n,u j ,u j _ 1} ...,Uf) + G(m,n,u k ,u k -i, ...u k >). Re- 
call that in the particular case OH]) this fact has been proved in [TJ] . Efficiency of the symmetry 
finding algorithm is demonstrated in §§4-5. In the fourth section a non-autonomous symmetry 

9 = (-l) m+n C 7 U{U2 ~ + 1} , (10) 

(UUI + U\ — U + l)(MM_i — U-\ + 11 + 1) 

is found for the equation 

u ltl u{ Ul - \){u x + 1) + (ui + l)(ui - 1) = (11) 



suggested in [SJ. In a recent paper [26] the connection of this equation with discrete SG equation 
is discussed. In section 5 a fifth order symmetry u t = g is evaluated for the same equation (fTTj) . 
where 



uu\ (u\ — 1) (u 2 — 1 



(ui + 1 + U\U - u) 2 (1 — Ui + u 2 + U 2 Ui) 



u~iu (u-i + 1) {u 2 - 1) ^ 

(— M-l +«-!«+ 1 + uf (1 + U-i - U_ 2 + W_ 2 «-l) 
U {u 2 — 1) (1 + MiM_i) (uU-i — U-i — UU\ — Mi) 



{U\ + 1 + UUi — u) 2 (1 — M_! + tt«_! + V? 2 

For some certain cases the systems fl2T|) . ([30]) . fl3Tl) . (jSJ) derived in section 3 are sufficient to 
find the final form of a fifth order symmetry if it does exist. It is the case for the equation ffTTT) 
mentioned above. However, they are not the only differential consequences of the integrability 
conditions listed in Proposition 1. In section 6 we briefly discuss how to derive additional 
systems of linear first order PDE consistency of which is also necessary for existence of higher 
symmetries. 

2 Higher symmetries and integrability conditions 

Existence of higher symmetries is believed to be an important indication of integrability. As- 
sume that chain ([1]) admits a higher symmetry (generally non-autonomous) of the form 

u t = g{m, n, Uj, Uj_ u Uf, u k , ...%), (13) 
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where j > 1, j' < — 1, k > 1, k! < — 1. As it is well known the function g should satisfy the 
linearization of chain (pQ) 

DDg = f Ul Dg + f Ql Dg + f u g, (14) 

Lemma 1. Linearized equation ( I14p can be rewritten in any of the following equivalent 
forms 

D-'Dg = fu^D-'g + f^Dg + f-^g, (15) 
DD-'g = ft-'Dg + ft'D-'g + f^g, 
D-'D-'g = f-^D-'g + f^D^a + fu 1 '-^- 

Proof. Let us prove the first statement of Lemma 1. To this end apply the operator D^ 1 to 
equation (fT4l and simplify it due to the relations 

D-\f u /f ai ) = -fZ\ 1 , D-\l/f Ul ) = f^\ D-\f ul /f Ul ) = 

The other formulas are verified in a similar way 

Differentiate equation f[T4"]) with respect to the highest order variable Uj+i and get (see [TTj ) 

DDg u .Dif Ul = f ul Dg Uj . (16) 

Set z : = log g Uj and rewrite the last equation in the form of a conservation law 

z 1 = z + {D- 1 -D^ l )\ogf Ul . (17) 
Differentiation of (TL5I) with respect to the lowest order variable uy-\ yields 

D '»<,„,D J 'f„ { y = f-^D 9u .,. (18) 
Obviously, equation (fT8|) implies for v := log^., 

^= u + (L>-^'+ 1 )log/-_ 1 ' 1 . (19) 

Actually, equations ( |T71) and ( |T9l provide necessary conditions of integrability of equa- 
tion ([1]): there exist numbers j and j' such that functions (D^ 1 — D^~ x ) log f Ul and (D — 
Di' +1 ) log/" 1 ' 1 belong to the image of the operator E — D, where E is the identity operator. 
These integrability conditions derived for the first time in [11] are highly nonlocal, since each 
equation contains the values of unknowns z, v taken at two different values of their arguments 
z = z(m } n,Uj,Uj_i, ...,Uj/ } Uk,Uk-i, ■■■Uk>), v = v(m, n, Uj, Uj-i, Uj>, Uk, Wfc-i, ■■■Uk>) and z = 
z(m,n+l,u j:1 , ...,Uj> t i,Uk+i,Uk, v = v(m,n+l,u j: i,Uj-i : i, Uf,i, Uk+i, «&, ...Mfe'+i). 

In other words equations (IT7|) and ffl~9l) are functional equations. 

To derive the next necessary condition of existence of a symmetry of the form ffT3l) with 
j > 2 differentiate equation with respect to the variable Uj 

DDig^jD^if^) + DD(g Uj ){Di(f u ) + D^f Ul )D^\f Ul )] = f Ul D(g u ._ 1 ) + (20) 

+f ai D(g u .)D^ 1 (f Ul ) + f u g u .. 

Put z^ 1 ' := g Uj _ x and rewrite equation (T2"Uj) as follows: 

zf' = r(m, n, U\, Uj, Uj-i, ...u-ijz^ + R(m, n, ux,Uj, Uj-±, ■■■Uf), (21) 
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where r = n-i^/i^l and 

R = Dj _^ ^ D- 1 {f ai D(g Uj )DJ- l (f Ul ) + f u9u . - DD(g Uj )Di(f u )} - D^D^f^). 

Thus a symmetry exists only when function R is in the image of the operator D — rE. 

Differentiating equation f|T5|) with respect to the variable uy one obtains for any j' < —2 

D-'Dig^jD^if-^) + D~ 1 D(g u ' j )[D j ' (f' 1 ' 1 ) + />''(/„ /' )/> r ' 1 i/„ V 1 ) = (22) 
= fu^D-'igu.J + fu^Digu^D^if^ 1 ) + ./•'••,/,,.. 

Put t?W := g Ujl+1 and rewrite the last equation as follows 

•yW = f(m, n, Mi, Mi, u, ...Uj' + i)v^ + R(m, n, Ui, Uj, Uj-i, ...ujr), (23) 
where the functions 



D(f-_\ 



and 

R = D fJ {f -Y) D {fn^D{g Uj ,)D3' +l {f-^) + I,, •'//,, - D ; /^//„,)/>M/, i 1 -)} - 

-D{g U]l )D^\f^) 

are expressed through known quantities f~ lfl , v, their derivatives and shifts. Obviously equa- 
tion (123]) provides one more integrability condition: function R should be in the image of the 
operator D — rE. 

Finalizing the reasonings above we obtain the following statement. 

Proposition 1. If an equation of the form ([1]) admits a higher symmetry of sufficiently 
great order then the conditions hold for some entire j and j': 

1) (D- 1 - LP- 1 ) \ogf Ul G Im(E - D) for j > 1 and f < -1 ; 

2) (D - Di' +1 ) log G Im{E - D) for j > 1 and j' < -1; 

3) R(m, n, Ui, Uj + i, Uj , ...Uji) G Im{rE — D) for j > 2 and j' < —2; 

4) i?(m, n, Mi, Uj+i, Mj, ...Uj') G Im(rE — D) for j > 2 and j' < —2. 

In such a way one can derive a sequence of integrability conditions given as functional 
equations. Obviously, in this form integrability conditions are not very effective. Below in the 
next section, we suggest a rule to derive some differential consequences of these equations which 
have a form of systems of first order partial differential equations and can effectively be studied 
by standard methods of differential algebra. 



3 Algorithm of finding symmetries 

Concentrate first on equation ( TT7|) . Note that z might depend only on the variables Uj, ...,Uj>, 
it does not depend on Uj with j ^ 0. Indeed, suppose that z depends on k > and does not 
depend on u s , s > k, then evidently z\ depends on Uk+i that is not possible since it contradicts 
equation (TT7|) . The case k < is studied analogously. Reduce (|T7|) to a system of the first order 
PDE's. To this end apply the operator D^ 1 -^ to both sides of the equation and have 

Y lZ = f<», (24) 
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where /« = D~ 1 -^(D~ 1 - L^ 1 ) log / Ul and Y x = D~ l -^D (see formula © above). Note 
that the coefficients of equation (|2^|) depend on the variable U-i, in spite of the solution z 
cannot depend on it. That is why one has to put an additional condition 

X x z = 0, (25) 

where Xi = -J?— . 

In a similar way the functional equation f lT7|) implies two more differential equations. Indeed, 
the equation can be represented as z_i = z — D^ 1 {D^ 1 — D^ 1 ) \ogf Ul . Applying the operator 
DX\ gives 

Y- X z = (26) 

where = -Y-x[{D~ x - D^ 1 ) log/ ttl ] and F_i = D-^D^ 1 (see formula © above). 
Thus we come up to a system of the first order non-homogeneous linear equations 

Xxz = 0, 

Y lZ = /« (27) 

X_!Z = 0, 

Note that for j = 1, f = — 1 system ( 127)) in essence coincides with that suggested in [TTJ. Com- 
patibility of system (127)) is necessary for existence of symmetry ( I13p . Emphasize that generally 
system ( 127)) is not closed. To close it we have to add all its linearly independent differential 
consequences obtained by taking cross applications of the operators such as [X 1? Yi]z = Xif^\ 
[Yx, Y_i]z = Yif(~^ — Y-if^ etc. In such a way we arrive at a system of the form 

L s z: = J2 a s,i^ = F {s \ s = l,2,...N, (28) 

i=j' 1 

X x z = 0, 
X^z = 0, 

satisfying the condition: any further cross application of the operators L s gives an equation of 
the form [L s ,L r ]z = L S (F^) — L r (F^ s ') linearly expressed through already known equations 
(125)) . Then the system is closed. 

Introduce notations = J^- and rewrite system ( 128)) as a system of linear algebraic 
equations 

Y f a B)i p® = F<- B \ s = l,2,...N. (29) 

i=f 

Two last equations in ( 128)) are valid automatically since z does not depend on u x ,U-x- Due 
to the well known Kronecker-Capelli theorem, system of linear equations (1291) is compatible 
if and only if the rank of the coefficient matrix A = (a S) i) is equal to that of the augmented 
matrix B obtained from A by adding the column of free terms F^ s \ Thus the condition 
rank(A) = rank(B) is necessary for existence of a symmetry f)T3T) . 

Let us deduce some differential consequences of equation (1T91 . Just applying the reasonings 
above to this equation we obtain 

X lV = 0, 

Yxv = p\ (30) 
X^v = 0, 
Y-xV = 
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where /« = D' 1 £.(D - Tfi log/" 1 ' 1 and = —Y_i[(D — D j ' +1 ) log/- 1 ; 1 ]. Similarly to 
the previous case one can close system (|30|) and then reduce it to a system of linear algebraic 
equations. We omit this part. 

Study the third condition of Proposition 1, which requires the consistency of the functional 
equation (T2"Tj) . Note that again unknown does not depend on the variable ui, hence applying 
the operator Z) -1 ^ to equation (|2~Tj) yields 

Yl z^ = D-\r ai )z { l\ + D~ 1 R Ul . (31) 
Exclude z_i by means of the same equation (12 ip to obtain 

Y x z<n = AWzM+B®, (32) 
where = D~ 1 (\ogr) ai and B^ = D~ l (R Ul — R(logr) Ul ). Rewrite equation f l2T|) as follows 

z^ = -(z?-R) (33) 
r 

and shift it back by applying the operator D~ x 

z { -l=D-\-)z^-D-\-). (34) 
Apply the operator Dq§-^ to the last equation to get 

Y^z® =Y- 1 (±)# ) (35) 
which obviously is rewritten as follows 

y_ x z« = A^z^+B^\ (36) 

where A^ 1 = rY_i(l/r) and 5 _1 = — Y_i(i2). As it was done above we put two more equations 
X\z^ = and X-\Z^ 1 ' = 0. Thus we find a consequence of the necessary condition 3) of the 
Proposition 1. The following system of the first order linear PDE should be consistent 

X lZ V = 0, 

Y_ lZ ^=A^z^+B^, (37) 
X_ lZ V = o. 

Deduce differential consequences from the fourth condition in Proposition 1, to this end 
apply the operator D -1 ^ to equation (|23|) : 

Y lV W = D-\r a Av { l\+D- l R ai . (38) 
Exclude v^l by means of equation (12"3"|) to obtain 

Y,v^ =A^v^ +B {1 \ (39) 
where A^ = D~ l (\ogr) Ul and B^ = D~ 1 (R ai — _R(logf) Sl ). Rewrite equation (|2"31 as follows 

v {l) = l(vP - R) (40) 



and shift it back by applying the operator D 1 

i£)= D-\\){vV - D~\^). (41) 



r 



Apply the operator Dq=—^ to the last equation to get 

r_ 1 ^ 1 ) = r_ 1 (^ 1) -r_ 1 (^), (42) 

which obviously is rewritten as follows 

Y. lV W =A^v^ + B(- 1 \ (43) 

where A' 1 = rK_i(l/r) and B' 1 = -lY^(R). As it was done above we put two more equations 
Xiv^ = and X_iV^ = 0. Thus we find a consequence of the necessary condition 4) of the 
Proposition 1. The following system of the first order linear partial differential equations should 
be compatible 

X lV {1) = 0, 

Y^ = A^v^ +B^ l \ (44) 
X_ lV {t) = 0. 

Continuing this way one can derive i-th integrability condition which is formulated as a 
compatibility conditions for a system of equations for the function := g Us with j' < s < j. 
It is easy to prove by induction that all of the functions g Ug depend only on the variables 
m, n, Uj, Uj-i, ...,Uf. This observation allows to prove the following statement on decomposition 
of a symmetry of the lattice (pQ): 

Lemma 2. Any symmetry ([TBI of the lattice ([TJ splits down into a sum of two functions 

g(m, n, Uj, Uj-i, iij>, u k , Uk-i, ■■■u k >) = F(m, n, Uj, Uj-i, %') + G(m, n, u k , u k -i, 

Here it is supposed that u = u = u. 
Proof. We observed above that 

g Us = a s (m, n, Uj, Uj-i, Uji), where f < s < j. 

These formulas yield 

9u s u r = where f < s < j, k' <r < k, r ^ 0, s ^ 0. 
The latter immediately proves the lemma. 



4 Evaluation of a non-autonomous symmetry for quad 
graph equation 

It is proved in [12J that equation (fTTl) does not admit any autonomous symmetry of the form 
(Q but it is conjectured there that it might have a non-autonomous symmetry: 

u t = g(m,n,u_i,u,ui). (45) 
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Apply the scheme above to equation (fTTj) to look for its symmetry (|45p . Function z = logg Ul 
solves equation ( !T7|) of the form Y\z = Y-\Z = p 1 '. Rewrite it in coordinates: 

u 2 i — 1 dz dz 2u\ dz u 2 — 2u — 1 

— — ^ 1 1 = ; 

2u dii-i du u 2 — 1 du\ u(u 2 — 1) 

2u-i dz dz u 2 — 1 dz u 2 ui + 2u 2 — U\ 



u 2 — 1 du_i du 2u du\ u(u 2 — 1) 

Note that the latter is not closed. To close this system one has to add one more equation 
[Yi,y_i]z = Yif^ — Y-ifW. Now solve the obtained system of three equations with respect 
to the partial derivatives of z: 

dz dz 2{u + 1) 



du-i ' dui ux(u + l) + l — u 
dz 2(«i-l) 1 1 1 

+ - + 7 + 



du ux(u + l) + l — u u u + 1 u — 1 

This triple of equations is an overdetermined system of equations for z which is consistent 
because the previous system is closed. Hence its general solution z is easily found. It contains 
arbitrary function Ci(m,n) depending on both discrete variables: 

, Cx(m, n)u(u 2 — 1) 
z = log ■ 



(uui + U\ — u + l) 2 

Substitution of z to ( |T71) yields log — 0- The last equation is easily solved Ci(m, n) = 

(— l) n C*2(m) where C^irn) is an arbitrary function of one discrete variable. 
Therefore by solving equation z = log g Ul we find: 

C 2 {m)u(u - 1) 

UUi + U\ — u + 1 

For the further specification consider v(u-i,u) = log(7„_ 1 = \ogg2, u ^ 1 which due to the 
general scheme solves the equations: 

u 2 _ 1 — 1 dv dv 2u 2 — 2uu-\ + U-x 
2u du-\ du u{u 2 — 1) 

2u_i dv dv 1 — u 2 — 2u 
u? — 1 du-x du u(u 2 — 1) 
Now solve them with respect to the derivatives 

dv 2(u - 1) 



du-\ u-\{u — 1) + u + 1 
dv 2(w_x + l) 1 1 

+ - + 7 + 



du u_i(u — 1) + 1 + u u u + 1 u—1 
and then find v : 

C 3 (m, n)u(u 2 — 1) 

v = lor 



(UU-1 — U-x + u + l) 2 
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Substitute it to the difference equation (TTTJ|) and obtain the relation 

-C 3 (m,n + 1) 
C 3 (m,n) 

solution to which is C^im, n) = ( — l) n C^(m), where C^(m) is an arbitrary function of m. As a 
result function g takes the form; 

g = -(-IT °^< U -V - ( _ 1} n C,(m)u(u + 1) + 

UU\ + U\ — U + 1 UU-i — U-i + u + 1 

By substituting the function obtained into the linearized equation and applying to the result 
the following operator: 

("-1 ~ l) 2 r -i 9 l D -id 1_ 

2 <9tii Mi <9wi u\ — 1 

we come to the relation 

1 1 7,2 _|_ 1 

- -g' 3 (u) + -ft,(tt) = — 2 (C 2 (m - 1) - C 4 (m)) 

which immediately implies: 

tfs'O) - -^ 3 («) + -*9M = 0, C 4 (m) = C 2 (m - 1). 
tt u 

Therefore g^(u) = Cs(m, n)u + Ce(m, n)u\ogu, where the coefficients do not depend on u but 
might depend on m,n. 
Now we have 

/ 1W C 2 (m)u(u-1) , . n C 2 (m - l)u(u + 1) , > , . 

# = -(-l) n f- - (-l) n — — fr + C 5 (m,n)u + C 6 (m,n)u\ogu 

Substitute it to the linearized equation and apply the operator 



Mi <9wi -uf — 1 ' 
to get: 

2(u 2 + l)C 5 (m, n) + 2((w 2 + 1) log w - u 2 + l)C 6 (m, n) + 2(w 2 - l)C 6 (m + 1, n)- 

-C 2 (m)(u - l) 2 + C 2 (m -l)(u+ l) 2 = 0. 

Collecting the coefficients before linearly independent functions (u 2 + 1) \ogu — u 2 + 1, u 2 — 1, 
2u, u 2 + 1 we obtain the system of equations: 

C 6 (m,n) = 0, C 6 (w + l,n) = 0, C 2 (m - 1) + C 2 (m) = 0, 

C 2 (m - 1) - C 2 (m) + 2C 5 (m, n) = 0. 

Thus 

C 2 ( m ) = -(-l) m C, C 6 (m,n) = 0, C 5 (m,n) = (-l) m (-l)"C. 
Due to the last equalities we find the final form of the symmetry desired: 

u(u 2 — l)(lliU-i + 1) 



g = (-l) m +"C 



(uui + ui — u + l)(nn_i — u_i + u + 1) 
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5 Evaluation of a fifth order symmetry for quad graph 
equation 

In this section we proceed with the testing of our symmetry finding algorithm. Apply it again 
to equation (II ip . Now we will look for a fifth order symmetry for it 

u t = g(u 2 , Mi, u, u-i, u-z). (46) 

Start with the equations (1241) and (126]) for the function z = \ogg U2 . These equations being as 
follows Y x z = /«, Y- X z = in an enlarged form are 

(n 2 x — l){u 2 _ 2 — 1) dz u 2 _ x — Idz dz 



4m_im du-2 2u du_i du 

2u\ dz 4m 2 mi dz 1 Au\ 



,2 



(47) 



w-l&U! {u\ - 1)(m 2 - 1) du 2 u {u\ - 1)(m 2 - 1)' 

4m_ 2 u_ x dz 2u_i dz dz 

(u 2 — l)(uL x — 1) du_ 2 u 2 — 1 du_i du 

uj-ldz (u§-l)(uf-l) dz 2u m 2 (m?-1) 

+ J—. IT" = 75 7 — • I 48 ) 



2u du\ Au\u duz u 2 — 1 2u\u 

Obviously system (|4"T|) . (1481) is not closed. To close it we add its differential consequences 
obtained by taking cross applications of the operators Y\, Y-\\ 

[Y^Y^z = YxfW - y_i/« =: /, (49) 
[Y 1 ,[Y 1 ,Y_ 1 ]]z = Y 1 f-[Y 1 ,Y_ 1 ]f^\ 

[y_ l5 [y^y^Wz = y_j - [y 1 ,y_ 1 ]/(- 1 ). 

System (|49p is closed and consistent. Solve it respect to the derivatives and find 

dz dz dz 2(ui + 1) 



dU-2 dU-l dU2 1 — Ui + U2 + U2U\ 



dz 1 2(1 - Ml ) 2m 
- + : + 



du u u\ + 1 + U\U — u u 2 — 1 ' 
dz 1 2(u + l) 2(1 -u 2 ) 2ui 



<9mi Mi Mi + 1 + MiM — M 1 — Mi + M 2 + M2M1 M 

Find z by integrating: 

MiM (m 2 - 1) (M 2 - 1) (M + 1) 

z = lor 



(1 - Mi + M 2 + M2""l) 2 (Mi + 1 + MiM - M) 2 ' 

which yields immediately an intermediate representation of g: 

MMi (Ui - 1) (M 2 - 1) 

= 7 — — ^27; ; ; r + 9i(ui, u, u_i, m_ 2 ) 

(Mi + 1 + M X M - U) (1 - Mi + M 2 + M 2 Mi) 
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For the further specification of g let us find v(ui, u, it_i, u_ 2 ) = logg u _ 2 = \oggi jU 2 . Func- 
tion v = v(u±, it, it_i, it— 2) should satisfy equations (compare with (I30I) ): 

(m 2 : - l)(w 2 2 - 1) dv u\ - 1 dv ^dv 



Au-iu du-2 2u du-i du 

M 2 - 1 9U! (« + 1) (U ~ 1) 2M_iU ' 1 J 

4w_ 2 M -i <9z 2m_i <9t> <9t> 



(u 2 — — 1) <9m_ 2 u 2 — 1 (9n_i (9m 

u\ — ldv 1 4u_i 



2u 9tti it + 1) (u-i -l)(u + 1) (u - 1) ' 



(51) 



System of equations ( )50l) . ( IBTl) can be closed by taking cross application of the operators 
Y\ and Y-\. After some manipulations of such a kind, which are omitted, one finds partial 
derivatives of v: 

dv 2(1 — w_i) dv 



<9it_ 2 1 + W-i — m_ 2 + U-2 u ~i ' dui 



0. 



<9v _ 1 2(w_ 2 + 1) 2(u - 1) 2m_i 



&U_1 W_i 1 + tt_i — M_ 2 + U_2 U ~1 —U_i + U_iU + 1 + It U 2 _ 1 — 1 

<9t> _ 1 2(u_x + 1) 2m 



(9m M — M_i + M_]M + 1 + W W 2 — 1 

Thus f is: 

M_]W ^M 2 ^ — l) (u 2 — 1) 



t> = log C 



(1 + W_i - M_ 2 + U-2U-l) 2 (-U-1 + + 1 + «) 2 



Since f = logg u _ 2 = loggi jU _ 2 by integrating we obtain more detailed representation for g: 

UU\ {u\ — 1) (w 2 — 1) 
g = ^ 1" 

(Ul + 1 + U\U —U) (1 — U\ + U 2 + U 2 Ui) 
_ M_iM + 1) (w 2 - 1) . . 

+^7 : 72— r + 92{u-i,u,ui), 

[-U-1 + M_iM + 1 + U) (1 + U_i - M_2 + U-2U-1) 

containing an unknown tail ^2(^-1, w, Mi). Due to equation (1321) function = dgijdu\ closely 
connected with z^ 1 ' = dgjdu\ solves the following system of equations 

u^-ldyW dyM 2u x dy^ (m 2 -1-2m) (1) _ 
2m du-\ du v? — 1 du\ u{v? — 1) ^ 

+ (m 2 -1)m 

(mi + 1 + MMi — m) 2 (— M_i + 1 + MM_i + w) 



2 1 



2m_! 9mW _ dyW + m 2 - 1 cfyW + 2 m 2 + mV - m 



v? — 1 du_i du 2u du\ u (u 2 — 1) 

M(Ml — l) 2 ^ 2 —2 M 2 M_i — M + 3 MM_i 2 + 2 M_i — 1 — M_ x 2 



U + 1 y (Mi + 1 + MMi — m) 3 (— M_i + 1 + tt«_i + w) 2 («i + 1 + MMi — u) 



12 



These equations imply 

0y« 4w 2 (w 2 - l) 2 (uiu_! + 1) 



<9w_l (— W-i + 1 + MM_! + m) 3 (Mi + 1 + UUi 



<9?/ (1) 2 (u + 1) (1) 2jn_i - 1) (m 2 - l) 2 -» 



<9«i Mi + 1 + MMi — It (u x + 1 + MMx — w) (—U-i + 1 + MM_! + «) 

We omit the expression for since it is huge. 

Thus we obtain the next representation for the searched g: 

UU\ (ui — 1) {u 2 — 1) 

g = 2 1" 

(ui + 1 + u\u - u) (1 - Mi + m 2 + U 2 Ui) 

+c M_in + 1) (m 2 - 1) ^ 



-M_! + + 1 + «) (1 + M_i - W_ 2 + W-2«-l) 

f— «i (u 2 — 1) ?i_! 2 + 2 (« — 1) (UUi — U + Mi) M_! + Ui (li + I) 2 ) (it — 1) U 

(—U-i + 1 + uu^i + u) 2 (ui + 1 + uui — u) 2 

C 2 - 1) % / v 

+ g 3 («_!,«), 



iti + 1 + ititi — It 

where C 2 is the constant of integration. By substituting this repsentation of g into linearized 
equation (1141) and differentiating it with respect to u 2 we get C 2 = 0. In order to look for g% 
differentiate the linearized equation respect to and study according to the scheme above 
the linear equations for = ^f 3 -- Their consistency condition immediately implies C = — 1. 
For this value of C the equations mentioned are of the form: 

u 2 _ x - 1 dy^ | dyW | w 2 w_i -2m 2 - u-i (2 ) _ (3 mm_i + u + w_i + 1) (m - 1) 



+ + — 7\ V 



2u <9w_i du (u 2 -l)u (w + l)(l — U-l + + m) 

2u_i <9?/ (2) <9?/ (2) + (w 2 -l + 2w) (2) _ - 1) (1 + u + u_ x - uu_i) u 



u 2 — 1 du-x du (u — 1) (u + 1) w (u + 1) (1 - n_! + mm_i + it) 

By solving these equations we find the next specification of g: 

UUi (ui — 1) (u 2 — 1) 
g = 2 1" 

[Ui + 1 + U\U —U) (1 — U\ + U 2 + U2«l) 
+ 1) (m 2 — 1) 



(— W-i + + 1 + m) 2 (1 + n_! - W_ 2 + U_ 2 U-i) 

U {u 2 — 1) (1 + MiM_i) — — MMi — Mi) (« + 1) wC 3 
2 2 1" 9i\ U )- 

(ui + 1 + mti - u) (1 - tt_i + + «) — w-i + 1 + + u 
By substituting it into linearized equation and differentiating the result with respect to the 
variable w_i we prove that C% = and that #4 (it) solves a linear homogeneous equation which 
has unique solution g± — 0. Thus the final form of the symmetry searched is (see also ffT2l) 
above) 

UUi (ui — 1) {u 2 — 1) 
u t = g = 9 h 

(■Ui + 1 + it^tt - u) (1 —U1+U2 + u 2 u 1 ) 

U_xU (■?!_! + 1) (u 2 — 1) 



'— U-l + U^iU + 1 + uf (1 + W-l - M-2 + M_ 2 M_i) 



tt 



(ti 2 — 1) (1 + «!«_!) — M_! — MoMx — U\) 



{u\ + 1 + uui — n) 2 (1 — m_i + + n) 2 
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6 Higher order characteristic vector fields and additional 
differential consequences of integrability conditions 

For some certain cases systems (1271) . (13"U1) . (1371) . (jSJ) derived in section 3 are sufficient to find 
the final form of the symmetry searched if it does exist. But for the other ones these systems 
do not provide enough information to find the symmetry desired or to make the final decision 
on its existence. Here we briefly discuss how to derive additional systems of linear first order 
PDE consistency of which is also necessary for existence of higher symmetries. To this end we 
can make of use the higher characteristic vector fields defined in such a way 

Y k = D~ k X_ 1 D k and Y_ k = D k X l D~ k for VA; > 1. (52) 

It follows from equation ( 1T7j) that z 2 = z + (1 + D){D^ — D J_1 ) log/ Ul . By applying the 
operator D~ 2 X_i to the last equation we find 

Y t z = f®, (53) 

where = D~ 2 -J^-(l + D)(D^ 1 — D^ 1 ) \ogf Ul . The coordinate representation of the vector 
field Y 2 is as follows 

Y 2 = X 1 + D-\Y 1 f)^ + D-\YJ- 1 ' 1 )^ + D-\YJ 1 )^ + D-\Y 1 fll' 1 )J- + .... (54) 
au\ ou-i ou 2 0U-2 

The coefficients of this equation depend on the extra variable U- 2 , hence we have to require 

X 2 z = 0, (55) 

where X 2 = 

Represent now equation (ITT]) as z^ 2 = z — (D^ 1 + D^ 2 )(D^ 1 — D^ 1 ) \ogf Ul . By applying 
the reasonings above one obtains a pair of equations 

Y_ 2 z = f ( ' 2 \ X^ 2 z = (56) 

where = -Y_ 2 (l + D^D' 1 - D^ 1 ) log/ UJ and X_ 2 = ^. 
Iterating this way one obtains a system of the form 

Y iZ = /«, (57) 
X t z = 0, —N <i<N, i=£ 0, 

which should be consistent for any value of the natural N. Here Xi = q§—, functions 
are expressed through the given functions /, f 1 '" 1 , f 1,1 and their derivatives and shifts. 
Recall that functions f^ 1 ^ are defined earlier in ( 1241) . ( 1261) . 

Obviously such kind extended consequences can be deduced from any integrability condition 
discussed in 82. 



7 Conclusions 

Integrability conditions for quad graph equations are studied based on the symmetry approach. 
Effective consequences of existence of higher symmetries are derived by using characteristic 
vector fields. These conditions can be successfully used for testing and for classification of 
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integrable quad graph equations. An algorithm of evaluating higher symmetries (generally non- 
autonomous) for autonomous quad graph equations is suggested. Efficiency of the algorithm is 
approved by applying to an example recently found in [9]. Fifth and third order symmetries for 
this equation are evaluated, the latter turns out to be non-autonomous. We have to emphasize 
that the main difficulty of non-autonomous case is the circumstance that all of the constants 
of integration, appearing when we solve the corresponding systems of PDE should now depend 
on the discrete arguments m, n. 
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